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Abstract 

The weak coupling limit for a quantum system, with discrete energy spectrum, coupled to a bose reservoir with 
the most general linear interaction is considered: under this limit we have a quantum noise processes substituting for 
the field. We obtain a limiting evolution unitary on the system and noise space which, when reduced to the sytem's 
degrees of freedom, provide the master and Langevin equations that are postulated on heuristic grounds by physicists. 
In addition we give a concrete application of our results by deriving the evolution of an atomic system interacting with 
the electrodynamic field without recourse to either rotating wave or dipole approximations. 



1 QUANTUM THEORY OF DAMPING 

Irreversible quantum evolutions now play a fundamental role in many areas of physics especially quantum optics. A large 
'"""'body of physical literature has been built up around the problem of describing in a stochastic model the effect of a source 
of quantum noise on a given quantum mechanical system, emphasising the quantum stochastic properties of the source 
I of the quantum noise. However, quite generally, the classical approach to the stochasticity which is taken ignores the 
^3 essential quantum nature of the problem. In this paper we wish to develop the weak coupling approach [1] to quantum 
2 damping. In a sense the noise fcilds used to model physical sources are still quantum in nature, as will be explained 
G later, and because they are arrive at from a well defined physical scaling limit do not require us to put in the desired 
features of the noise by hand. 

The derivations of the quantum master equations and quantum Langevin equations for open quantum systems which are 

present in the current physics literature are well motivated from the physical point of view, cf. [2] , however mathematically 
j-- — imprecise. The heuristic procedure to render the reservoir, to which the system is coupled, into a source of quantum 
^sO noise via some markovian approximation is generally ill-defined and reliant on second order perturbation theory. 

On the other hand, the weak coupling limit for an open quantum system gives a device for obtaining irreversible evolutions. 
O^l Mathematically rigorous derivations of the Langevin and master equations along these lines have been given for certain 

Specific models by Pule [3] and Davies [4] . Other attempts have been made to approximate a quantum reservoir by fitting 
y—{ Quantum Brownian Motions (QBMs), cf. [5], [6] and [7]. 

^ In a long series of papers [[1], [8], [9], [10], [IT]], Accardi, Frigerio and Lu have developed their approach whereby a quantum 
reservoir can be reduced to a quantum stochastic noise source via a weak coupling limit procedure. The theory is 



. mathematically rigorous while at the same time applicable to the wide range of phenomena considered by physicists and 



C3 gives a precise description of the reservoir as a quantum noise source. The main mathematical device for establishing 
convergence of the reservoir feilds is a quantum central limit theorem; that is a central limit theorem for non-commuting 
observables: the theory of non-commutative probability, or Quantum Probability as it is more correctly known, affords the 
necessary mathematical framework to interpret the limit process. Indeed the Fock space descriptions of a bose reservoir 
and a quantum stochastic processes are so similar in nature, it is not surprising, with hindsight, that the former may be 
reduced to the latter following some limiting procedure. 

Our objective in this paper is to review the results of the programme of Accardi, Frigerio and Lu so far and to extend 
their theory so as to deal with the most general interaction between a system and a noise source encountered in physical 
theories. This we do and show that the energy shifts, linewidths, master equations and Langevin equations, arising for 
the system as a result of its coupling to the noise field, concur with those obtain by earlier researchers [2]. However we 
have, in addition, also a quantum stochastic description of the reservoir noise fields themselves. As a concrete application 
of our theory we consider the particular case where a quantum clectrodynamical field acts as reservoir, however we stress 
that this is only one of the many applications of our theory. 
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Wc shall discuss only minimal coupling intcracions, that is interactions linear in the creation/ annihilation operators for 
the reservoir. In a later paper we discuss how to treat the situations where the interaction is of polynomial type. 



1.1 Open Quantum Systems 

Wc consider a system (S) coupled to a reservoir (R). The system (S) is to be a quantum mechanical: its state space 
shall be a separable Hilbert space Jig. The reservoir, on the other hand, is comprised of one or indeed several quantum 
fields, and so has infinitely many degrees of freedom. We shall consider a bosonic reservoir; the state space for (i?) is 
the bosonic Fock space T-Lr over a separable Hilbert space H R ; in standard notation we write TLr = TB{T~i R ). H R is 
again to be a separable Hilbert space and may quite generally describe not only one but several individual species of 
particle in the reservoir. For instance, consider several species of particle Pi, P2, P3, ... in the reservoir and suppose that 
H R = ®j~Hp., where "Hp. is the state space for particle type Pj, then 

Ur = T B (H R ) = T B {®jH} Pj ) = ® 3 T B (H Pj ). (1) 

The space Tip is referred to as the (combined) one particle state space for the reservoir. The overall state space for the 
combined system and reservoir is Hs (DHr- The vacuum vector of the reservoir space will be denoted throughout as 'f p. 
In the following we consider only bosonic species in the reservoir, however it is also possible to work with fermions [5]. 
The dynamics of the combined system and reservoir is governed by the formal Hamiltonian which we may write as 

H W = H (o) + XHl . ( 2 ) 

that is, as the sum of a free Hamiltonian H^ and an interaction Hi, with A a real coupling parameter. is to be 

expressible as 

= H s ® 1r + l s ® H R , (3) 

where Hs and Hr are self-adjoint operators on the spaces Us and Hr respectively. For each A, we consider the unitary 
operator on % s ® defined by 

V t W =exp^H^. (4) 

This gives the time evolution under H^. A standard device in pertubation theory is to transform to the interaction 
picture; this involves introducing the operator 

Df^e (5) 

is a unitary operator on Hs ® T~Lr called the wave operator at time t. We note that {U^ : t e 1} is a left 
-cocycle, that is it satisfies the relation 

U$ s =v?\u^)U^. (6) 
We define the time-evolutes for X 6 B(Hs) <8> B(Hr) as 

V M {X) = V W XV M, 



u 



W {X) = U ( t X)] XU < t X) . (7) 



The Schrodinger equation for the time evolutions are 



and 



1 V W = L H W V W 
dt * ih 1 



With these we deduce the associated Heisenberg equations 

d 



-V<?\X) = h4 X \x)M X) (HW)] = U X \[X,HW}) 
at in in 

and 



§ i u?\x) = ±ul»([X,v?\H I )]). (9) 



From (8) we obtain the integral equation 



ih 

and consequently the iterated Born series 

00 \ n nt rtn-i 



C/ t (A) = l + 4 fdsvpHHriuW, (10) 







1 + E(^ / . tfl-/ dtn^CHj)...^^/). (11) 
n=l * ^° *' 

This may be resummed as 

[7-W=Texp{A ^dsv^iHj)}. (12) 
«fi Jo 

where T denotes time ordering. 

Before we continue, we must say more about how to interpret the formal sum of and Hi. Firstly we assume that 
and Hi are self-adjoint operators on Hs <8> Wfi- We shall assume that, for sufficiently small A and bounded i, the 
iterated series (11) is uniformly convergent and is bounded on the algebraic tensor product of Us nd the 

set of exponential vectors. From the cocycle relation (1.1:6) we have that if we define the unitary operator by 

V t W = V t {0) lll X) (13) 

then {V t W : f £ 1} gives a strongly continuous unitary group whose generator j^H 1 -^ is formally given as 
H(°> + XHj. The time evolution in the Heisenberg picture is then given by 

v?\x)=vi X \v<?>(X)). (14) 

1.2 The Free Evolution 

For U\ = L 2 (K 3 ) (momentum space), we have creation/ annihilation operators a"(k) satisfying 

[a(k),a(k') + ] = <S(k-k'), [o(k) t ,a t (k')] =0, [a(k), a(k')] = 0. (15) 
The creation/ annihilation fields on are defined, for G H R , by 

At( 5 ) = /" d 3 fcff(k)a f (k), A(ff) = /" d 3 fcp(k)a(k). (16) 

From (16) we obtain the canonical commutation relations (CCR): 

[A(h), A'(f)] = (h, /}, while [A(h), A(f)] = 0. We take Hr to be the second quantization of an operator H R on H R given 

by 

(Hym = hw{k)f{k). (17) 

H R may then be expressed as 

H R = / dkHuj{k)a f (k)a(k). (18) 



Note that v t (0) (l s ® A*{g)) = 1<? <g> A*(S t g); whe re we have introduced the unitary operator S t on % R given by 



S t =exp{~H R }; (19) 



that is (S t f)(k) = eM k ^f(k). 



1.3 The standard approach to the quantum Langevin equation 

Consider an interaction of the type 

H I = ih{D®A\g)-D^®A{g)}. (20) 

where D € B(l-Ls) has a harmonic free-evolution; that is — \D,Hs] = —iuiD. Then we have 

in 

vi \Hj) = ih{D ® A\S?g) -X>t ® A(S?g)}, (21) 

where 

S? = e~ luJt S t , (22) 
that is,(S£7)(k) = e^^-^'/O). Now write X t = u[ X) (X <8 1 R ), for X E B(Us), then from (9) we have 

^1 = \{A\{S?g)[X,D] t - [X,D%A t (Sfg)}, (23) 

where A\ (/) = ^ A) (ls 8 ^(/))- 

Now one can show that A t (f) = A J °° dsD s {f, S^g) + Is (8 ^4(/), so substituting back into (23) gives 



9X 4 _ 2 



/ {Dt^(s - t)[X, D] t - [X, D%^(t - *)£>.} + \{&[X, D] t - [X, D%{f}, (24) 
Jo 



where 4>"(t) = (g, Sf_ t g) = J d 3 fe|p(k)| 2 e _i ^ (k)_a,)t and = l s ® In standard terminology 0(f) is called the 

memory function and the fluctuating quantum force [2,17] or input field [18], albeit in the interaction picture. One 
notes that, in the vacuum state, is gaussian distributed and all first and second moments vanish except the two-point 
function 

(*«,erer + *fl) = (s?g,s? g ) = <i> u (t-*)- (25) 

The standard approach taken at this juncture is to introduce the so-called first Markov approximation. Here, for example, 
one takes 'H R = L 2 (K), g = \ — (constant) and u>(k) = k. Then 

</> u {t) = / dk—e l{k+ul)t = Ke iut 5(t). (26) 

J-00 27T 

There are, however, several important objections to be made to this approach. Firstly any physical details specific to the 
reservoir must be put in by hand. Secondly, the spectrum of H R is here unbounded below, this is necessary to produce 
the delta function correlation of white noise. From a physical point of view this is unacceptable as H R must be bounded 
below for stability. Finally, the fact that the frequency spectrum ui(k) = k is unbounded below precludes any possibility 
of dropping the rotating wave approximation. 

1.4 The Weak Coupling Limit 

We now describe the ideas behind the weak coupling limit in the simplest situation where we have taken a dipole and 
rotating wave approximation. We define the following reservoir operator 

B { t ^ x \g) = A(X f dt x S£g) =\f ^{dk) f dt ie -^-^g{k)a{k). (27) 
Jo Jm Jo 

Calculating the two-point vacuum expectations gives 

(* R (0),Bt x \g)B^(m R (0)) = \ 2 fdh f d Sl (S-g,S-J) 

Jo Jo 

du [ dT(S?g,f), (28) 

Ju/X 2 -s 

where we have substituted u = \ 2 ti and t = t\ — s\. 



In order to obtain a non-trivial two-point function in the limit A — > we must rescale time as 



t^t/X 2 ; (29) 
this is known as the van Hove or weak coupling limit in physics. One finds 

/oo 
dr(S^gJ). (30) 

Physically the limit A — > with t ^ t/\ 2 allows us to consider progressively weaker interactions which are allowed to 
run over increasingly larger periods of time and so we obtain the long term cumulative effect of the interaction on the 
system. Now the creation and annihilation operators are gaussian in the vaccum state and, as a result, so too are the 
operators B^J^'K 

Furthermore the limiting two-point function (29) is suggestive of the correlation function of a Brownian motion. However, 
an interpretation of the above in terms of classical Brownian motion is erroneous as it ignores the essentially quantum 
probabilistic nature of these processes. 

1.5 The Interaction 

For technical reasons we work with a system Hamiltonian Hg which has discrete spectrum. is taken to be bounded 
below as required from physics. The type of interaction Hi which we wish to study is of the form 

H i = in E E i D 1 ® A ^9l) ~ Df ® A(tf) }, (31) 



where J 7, is a discrete subset of R. For each u> S F, we take 6 B(TLs) to have harmonic free evolution with frequency 

-q.[D*,Hs] = -ujD*, j = l,...,N(u). (32) 

Thus the superscript ui labels harmonic frequency and j = 1, ...,N(lj) the degeneracy of that frequency. An interaction 
similar to (30) has been treated in [10], however there the test functions gf were taken to be equal for each value of lo. 
Our reasons for studying (30) above are because it allows us to treat the most general interactions encountered in physics. 
Typically in quantum field theory one considers an interaction of the type 

ffj = ih J d 3 k {0(k) ® a f (k) - eHk) ® o(k)}, (33) 

where {^(k) : k £ R 3 } is a family of operators on TLs- The operators 9(k) are called the response terms: they contain 
local information about the interaction. In the dipole approximation of quantum field theory one makes the replacement 

0(k) dipolc (k) = g(k)6(0), (34) 

where g(k) is some suitable test function. The physical argument is, cf. [2], that the response does not vary appreciably 
for values of the wavelength of the reservoir particles which are large relative to the physical dimensions of the system, 
though this can hardly be true for large momenta. As a result one obtains the approximate Hamiltonian 

H z Hf polc = ih{9(0) ® A\g) - 0\0) ® A(g)}. (35) 

A further approximation often made by physicists is to replace 9(0) by an operator D having a harmonic free-evolution 
with some frequency uj £ R. This approximation is just the rotating wave approximation. 

In order to avoid these appproximations we argue as follows: Let B be a complete basis of eigenstates of Hs, then 



but however we may write 



Hj= ]T {4>\Hi\(j)')\4>){(t)'\, (36) 

4>,<t>'£B 

(<t>\HiW) = ihl dkmewwWi*) - (0i^(k)i0V(k)} 



= in[A^^)-A{g^)], (37) 

where we have introduced the test functions 

ff ^Ck) = (0|fl(k)|^): (38) 

note that the order of (j> and <f> is reversed in the second term in due to the conjugate linear nature of the creation fcild. 
This now means that the interaction can be expressed as 

Hj = ih J2 { T W® A ^9W)- T U>® A (9W)}, (39) 

<p,<f>'£B 

where we have introduced the transition operators Ta^i = |<^)((^'|. 

We note that the transition operators T^i are harmonic under the free evolution, in fact we have 

tt , H s ] = -iu^y T u > , (40) 

where 

Em — E<f, 

= h ■ ( 41 ) 
So F = {ui^i : 4>, 4> S B} is now the set of Bohr frequencies. 

The expression (38) is now equivalent to the interaction (30) which we propose to study. Here we need only relabel the 
as where uj = uj^i and the j again labels degeneracy. The functions are relabeled accordingly. 

2 THE QUANTUM STOCHASTIC LIMIT 

2.1 Quantum Brownian Motions 

In this section we first of all discuss the concept of quantum brownian motion. As this is not yet widely known amongst 
physicists we give an exposition below: 

Definition. A quantum brownian motion (QBM) is a triple (H,^>, (B t )t), where H is a separable Hilbert space, $ G H 
with ||$|| = 1, (Bt)t is a family of operators on H such that 

(i) qt = ReBt and pt = ImBt are classical brownian motions for the state 

(ii) [p s ,qt\ = — min(s,t), where k G R. 

The basic example is the following: Let He = Tb(L 2 (K)) and $c be the vacuum state. Then define B t to be 

Bt = A c ( Xm ) (42) 
where Ac is the annihilation operator on He- From the (CCR) we have 

[B U B\] = (x [0)t] ,X [0iS] ) = min(M); [B t ,B a ] = = [flJ.Bj]. (43) 

So setting q t = —(B t + B\) and p t = — (B t — b!) wc have from the (CCR) that [pt,q s ] = —[itiPs] an( l 
2 2i 

\pt,q s ] = imin(i,s). (44) 

Now if wc set 

dBt = B t+dt -Bt = Ac(X[t,t+dt]), (45) 

we have that 

($ c ,dB|$c> = 0; ($c, (dBl)H c ) = 0; ($ c ,dBldB t $c) = 0; 

($ c ,dB t dBl$c) = dt. (46) 

Now dB t and dB\ are gaussian in the vacuum state, because the creation and annihilation fields are, therefore so too are 
dq t and dp t - Furthermore, 

<$c, (dq t )H c ) = i($ c , (dB t + dBt) 2 * c ) = \dt (47) 



and similarly ($c> {dpt) 2 ^c) = ~^dt. Finally noting that at unequal times s and t 

(<S>c,dq t dq s <P c ) = = ($ c ,dptdp s $c), (48) 

whenever t(t + dt < s(s + ds or s(s + ds < t(t + dt, we conclude that (qt)t and (pt)t are each separate Brownian motions 
for expectations taken in the state $c- So {He, &c, {Bt)t} is a QBM. We can introduce formal creation and annihilation 
densities 6"(t) satisfying 

[&(*), 6(*)] = = [6 + (t), [6(t), &+(*)] = S(t - a), (49) 



such that 

At 



(g) = [ dsg( S )b( S ); A^g) = f ds g(s)b + (s); (50) 

JR JR 

from this we see 

B\ = [ dsb $ (s). (51) 



dB 

We may write b\ = — - and consider these densities as "quantum white noises" . 
dt 

More generally let K be a separable Hilbert space and let L 2 (R, K) denote the set of square-integrablc K- valued functions 
over R. Now h £ L 2 (R,K) is a function h(t) G K with J R dt \\h(t)\\ 2 K (oo. The inner product on L 2 (R, K) is given by 

(h,ti) = [ (h(t),h'(t))dt. (52) 



If {e n } n is a complete orthonormal basis for K then we can write h(t) = h n (t)e n , where h„(t) = (e„, /i (£))#-; this 
gives a natural isomorphism 

L 2 (R, K) = K ® L 2 (R). (53) 

Now take TLk = rs(L 2 (R, K)) and let <£># denote vacuum vector of Hk, then a QBM is given by (Wki&k, (B t (g))t), 
for non-zero g 6 K , where 

Bt{g) = A K {g®X[o,t)) ( 54 ) 
where Ak is the annihilation operator on T-Lk- The commutation relations are 

[B t {g),B\{f)] = (s®X[o,t]./®X[o,»]> = (ff./>JC ™ii(t,«), (55) 

with remaining commutators vanishing. So {H_r- = rs(L 2 (R, -?0), ^if , (Bt(g))t} i s a QBM. Taking K = C and |g| 2 = re 
leads back to the original example. 

However, there is a more general possibility than that above; Let Q > Ik and set 

C = Q®1 L 2 (R) . (56) 

Then let ip c be the state on the Weyl algebra W(TLk) with covariance C. We can construct {Gb(Hk,C), 7t^ 1k ,^^ k } 
the GNS triple over {W(Hk), C} and on it define the operator 

B Q (9,t)=n° K B t (g). (57) 

Then {Gb(Hk, C), $>^L > {Bq{9, t))t} is a QBM referred to as quantum Brownian motion over L 2 (R, K) with covariance 
C, or more loosely with covariance Q. We have that 

^07,*)Bj ) (/,*)*S Jf ) = M5 t 07)5i(/)) =mm(t,s)(g,Q±±f) K , (58) 

and similarly 

($£ K , B* Q (f,s)B Q (g,t)*% K ) =min(t,s)(g,Q—±f) K . (59) 



2.2 Quantum Stochastic Calculus 

As is well known, a stochastic calculus can be built up around classical brownian motion and that the resulting theory has 
widespread applications to studying noisy systems in physics and engineering. It is also possible to build up a quantum 
stochastic calculus based on the QBMs we have just considered. This was originally done by Hudson and Parthasarathy 
[T5] . [14] , The basic integrators are dt and, depending on the context, dB\ or dB\(g) or dBg(t,g). 
In the simplest case, for instance, we have, for a partition — oo = ti < < 2 < ... < t n < t n+ i = oo, 

n 

L 2 (K) = 0L 2 ([i m ,t m+1 ]) 

m— 1 

and consequently 

n 

r B (L 2 (K)) = (g) T B (L 2 ([t m , t m+1 \). (60) 

This gives the required time filtration in the quantum situation. We say that a family of operators (X t ) t on T b {{L 2 (M.))) 
is adapted if, for all t, 

X t = X t ®l (61) 
on Tb(L 2 ((— oo, t)) <S> rs(L 2 ([i, oo))). The quantum Ito table reads as 

dB Q (g,t)dB^ Q (f,t) = (g,Q±±f) Kdt; 

dB f Q (f,t)dB Q (g,t) = (g,9^1f) Kdt} 

(dt) 2 , dtdB* Q (g,t), (dB* Q (g,t)) 2 = 0, (62) 
Let (X t )t be an adapted process and of the form 



X, 



J (x s ds + xfdB^g, s) + x~dB Q (g, s)) (63) 



and (Y t )t a similar process, then we have the quantum Ito formula 

d(X t .Y t ) = dX t .Y t + X t .dY t + dX t .dY t (64) 

with 

dX t = x t dt + x+dB^ Q (g,t) + XtdB Q (g,t)). (65) 
2.3 The Weak Coupling Limit of Quantum Field Theory 

The results of Accardi, Frigerio and Lu concerning the weak coupling limit for an interaction (31) which has undergone 
both a dipole and rotating wave approximation can be summarised as follows; 
Recall that 

/oo 
dT(S?g,f); 

now define a scsquilincar form (.|.)" on 'Hr R the one particle reservoir space by 



(g\fr= / dT(S?g,f). (66) 
We consider the space of suitable test-functions T w C H^, determined by the condition 

dt\(g,S?f)\(oo, (67) 



whenever f,g G T u . Note that technically T u does not depend on w however we keep it in as a label. Then we construct 
K w the completion of T u with respect to (.|.)". That is K u is the completion of T" factored out by its (.|.)"-norm null 
space. K 1 ^ is a separable Hilbert space with inner product 



Theorem 1. In the limit A — > the stochastic process on the resevoir space 

(a,, 
t/> 



for f € K u , converges weakly in the sense of matrix elements to QBM on L 2 (R, K u ). We denote this QBM by {H u = 
Tb{L 2 {R 1 K u )),G> u = (Bt (/))t} In the next theorem we show that uty 2 converges to a stochastic process U t on 

Hs <8> H u in a sense to be made explicit now. 

Theorem 2. Let /«, ftO") e A'"; T^, SC?")}0, for j = 1, n : j' = 1, m and Zet 0, </>' G i? s then the limit as A -> 
o/ i/ie matrix element 

<^4u)^(/ (1) )-"^ (68) 

exists and equals 

(0 ® (/«)...B#„, (/(«))$"|[/ t |0' ® B^ (ft(D)...B^ m) (ftH)$"), (69) 
where Ut is a process on Hs ® -ff" which is the solution to the quantum stochastic differential equation 

dU t = {D® dBf{g) - £>+ ® dB t w ( ff ) - (ff|ff) w -D+D <g> dt}tf t) (70) 

(g\fT-=[ d T (g,S?f). (71) 



Note that d(UtUj) = = d{U}Ut) by the quantum Ito formula and the Ito table. So Ut is unitary on ^jgM", however it 
describes an irreversible evolution when restricted to Hs- The unitarity condition corresponds to a fluctuation-dissipation 
law. 

Theorem 3 Let X 6 B(Hs), then in the notation of theorem 2 the limit 



Ul%? (X ® l fl )[$> a |0' ® fl^J a (h^)...B^J x , (72) 

exists and equals 

<^ (8) (/W)...^ 3 (/^^^ICT-tCJT ® l)l7 t |0' ® SQC^ 1 ))...^^™))*"). (73) 

Note that in these theorems we encounter vectors of the type B^,'^ (f)^ R which are exponential vectors with test 

functions A drS^ f: they are called collective coherent vectors in the terminology of Accardi and Lu and they are 

designed to extract the long time cumulative behaviour of the reservoir fields. 

2.4 Non-Zero Temperature Reservoir 

Next, for the non- vacuum case, we consider a density matrix pg on H-r which is invariant under the free evolution and 
gaussian with covariancc Q > l^i . That is 

1 

Tr{p Q W(g)} = e -2 <9 ' Q9> , y g e U\. (74) 

The invariance condition is equivalent to 

[S t ,Q] =0,onDom(Q). (75) 
In particular, the choice of a heat bath at inverse temperature (3 and fugacity z is given by 

1 + ze-P H R 

^ ~ i -8H 1 

1 — ze p r 

that is 



(76) 



(Qf)(k)=coth^(fkj(k)- f i)f(k), (77) 



where /i = — lnz is the chemical potential. 

P 

Now 

/oo 
JrTr{ PQ A(S-g)AHf)} 

f°° O + 1 

= min(M)y dr(5- 5 ,(^-)/), (78) 

and similarly 

Urn Tr{pg B^t^jB^AJ^j = min(t) s) £° dr( <^ ; (^)/>. (79) 
Let Tg be the subset of Dom(Q) such that 

/oo roo 
\(f,S?h)\dt <oo<md \(f, S?Qh)\dt < oo, (80) 
-oo J —oo 

whenever f,hE Tq. Let i^g be the Hilbert space completion of Xg with respect to the sesquilinear form (.|.)q given by 

/oo 
(f,QS?h)dt. (81) 
-oo 

Note that in most cases Tq is dense in and that i^g is a Hilbert space equipped with inner product (.|.)q- 
Theorem la: The process (i?|^ 2 (/)) t in £/ie mixed state Pq converges weakly in the sense of matrix elements to QBM 
over L 2 (R,Kq\) with covamance Q. This is denoted as {H^ = G B {{L 2 {R,K%),Q ® 1), (B%(f,t)) t }. 

Theorem 2a; For <t>,<j>' £ H s , fV\hV"> g Jf« and T^,S^)0, for j = 1, ...n; j' = l,...,m, £/ie fimii osA^Oo/ 

Tr{ \<j>')(<j>\ ®5^/ A2 (/^ (82) 

exists and equals 

(4> ® B2 t (/ (1) ' T(1) )- 5 O t (/ (n) ' T(n) )*Sl c; *l^ ® B$(h< 1 \sW)...B$(h<' m \sW)&Z) (83) 
where Ut is a unitary operator on Hs <S> G b{L 2 (R, Kq), C), with C = Q <g> 1, satisfying the QSDE 

dU t = {D ® dSg^o, t) - ® dBgfo, t) - (g|g)g+ £> f D ® dt - (fl|fl)JlDD t (8 di}*7 t , (84) 

w/iere 

(ffl/)g± = fj T ^^r^f)- (85) 
Theorem 3a: For X £ B{Hs), then in the notation of theorem 2 

hmJr{|0O<0l®4<^ ( 8fi ) 
exists and equals 

(<l>®B$(jM,TM)...B$(jW,TW)^^ (87) 



2.5 The Quantum Stochastic Limit for the Full Interaction 



Now suppose that the interaction is of the form (31). The problem of dropping the rotating wave approximation was 
first tackled by Accardi and Lu in [10] for an interaction similar to (31) except that all the test functions were taken to 
be the same. The result is that for each Bohr frequency u> we obtain a separate independent QBM. 

First of all note that r B (0 wgF L 2 (R, K%)) = <g> weF T B (L 2 (R, K%)) then consider the Weyl algebra W(0 weF L 2 (R, K%)) = 
<SWf W(L 2 (R, Kg)), with quasi-free state ipg with covariance C = (g) weF C, where C = Q ® 1 on U K % = K% ® L 2 (R). 
The GNS triple over {W(® weF L 2 {R, K%)), (p d } is = T B (® weF L 2 {R, K%)), <?)), 7rg, Now observe that 

Hq = «Q> «Q = (g) ®Q = (g) *q- (88) 
uef uef mef 

For each / G Kg we have 

BZ(J,t)=*%A%(J®X[o,t]), (89) 
where Ag is annihilation operator on L 2 (R, -Kg), so for t w )0, / w G -Kg, for each w e f , we have 

^((g)/ W) (U-eF) = (g)^(/ W) ^). (90) 



Theorem lb: For each uj G F and f G Kg the limit A — > 0. (/) taken in the state pg converges in the sense 

of matrix elements to QBM over L 2 {R,Kg) with covariance Q and each of these limiting processes are independent for 
different values of u>. 

Theorem 2b: Let fi 3 \h { P G K%;T$\S^)Q for each wGfj = 1, ...,n u ; j' = l,...,m u ; and t > 0; (f>, 4> G H s then 
the limit as A — > o/ i/ie matrix element 

Tr{|0')(0| ® [(g) ^^ /A2 (^ ) )---^L(^ 1) )] t x ^ a2 (/W)...^, a2 (/(""))] E$[ a } (91) 

exists and equals 

('!>»[^B^aL 1 KTP)...B^(fM,T^)]^\ (92) 
®[®^(/»£ 1 ^S£ 1) )...^(^^^)]t$g>, (93) 
where Ut is unitary on Hs <S> Hg and satisfies quantum stochastic differential equation 

du > = i E E W ® - *) - D f ® (-9" . *)] 

Ui£F j=l 

-EE («)«+*-£ E W f WI#5£ *> tf*. ( 94 ) 

wif/i f/o = 1. TTie coefficients are given by 

(f\hr Q ±= f dt{f,S?(^)h), (95) 

/or f,h<E K u . The Ito table is given by 

dB%(f,t)dB$*(g t t) = 5 UtU ,(f\g)% + dt 

dB^(g,t)dB^(f,t) ee 8 UiU ,(f\g)%_dt. (96) 

The coefficients in (2.5:8) are 

,Q + 1 



(f\g)Q+ = / dt(f,sn^)g) = (f\g)Q- + + W)^ 



(f\g)Q- = r dt(f,sn^)g) = (g\fr Q z + (f\g)%i. (97) 

The proof of theorem 2a is as follows; first of all we know that for different to we can always consider independent 
Q-quantum Brownian motions. This is done in [10]. The next step is to consider the effect of the degeneracy which 
may arise for each lo £ F. In this case we must, therefore, generalize the results of [8], [10] accordingly. This is done in 
appendices B and C. 

Theorem 3a. In the notations of theorem 2a, for any X £ B(Hg), the limit as A — > 

W>M ® [(g) ^a^ 2 (^ 1) )---^ ) )/A2 (^ ) )] t P Q [<g) B ^ /x M ) )-^J )/a2 (/^ ) )] u§l a (X ® 1)} (98) 

exists and equals 

(0® [(g) B^(fi 1 \T^)...B^(f^\T^^\uHx ® 1)W ® [(g) i?^(^ 1) ,5( 1 ))...^(^),^))]tcf>g), (99) 

lj£F ujEF 

where Ut is the solution to the quantum stochastic differential equation of theorem. 



2.6 The Langevin and Master Equations 

In each of the cases the right hand side of the expression for dUt contains a term of the form —Y ® dtUt- For instance, 
we have 

y = J2 E< D f D % W)q+ + D " D k WWT Q _ h (ioo) 

u£F i,k=l 

The Langevin equation then reads as follows 

d[u}{X®l)Ut] = [dU t y(X®l)U t + ul{X®l)dU t + [dU t }' [ (X®l)dUt 

= Ut [L (X)®dt+J2 E L 1+ W ® dB Q C# , t) 

ueF j,k=i 

+ E E L U X ) ® dBgOtf.*)}^. (101) 

wGF j,fc=l 



where 
with 



£ (x) = -xy - y*x + e(x), (102) 



N(oj) 

e w = E E{^ + ^[«i«+ +WsW Q ^ +D l XD kM^ Q - + wr Q --]} 

uieF j,k=l 

N(ui) 

= E E { Sf*^ + D%XDf {g^\gtr Q .\ }, (103) 

u)£F j,k=l 

and 

£?+(X) = XDJ - D^X, Z£_(X) = D^X - ID 3 wt . (104) 

Note that unitarity follows from 



Lo(l 3 ) = ~(Y + yt) + 9(l s ) = 0; £^±(ls) - 0. (105) 

1 i 

It is instructive to set Y = —T + j^H' s where both T and H' s are self-adjoint; we then have that 

l (x) = -hxr + rx) + e(x) + i[x, h' s }. (ioe) 

2 



The unitarity condition is then 2Rc Y = T = 6 (Is); this is the fluctuation-dissipation relation. The presence of the 
imaginary term H' s does not effect the unitarity. For pg a, density matrix on we define the expectation (.) t by 

(X) t = Tv{ Ps ® U}(X ® l)U t ] = Tr{ St X}, (107) 

where s t denotes the effective density matrix on (S) and the second trace is a partial trace over the system space (the 
trace over the reservoir space assumed to be taken already) . Now 

j f (X) t = Tr{p s ® jU}{X ® l)U t }, (108) 

so in terms of the effective density matrix St 

^Tr{ St X} - Tr{s t L (X)} = Tr{L* (s t )X}, (109) 

where Lq denotes the adjoint operation to Lq on the dual of £>(%s). This gives the master equation 

^ = L*(st) -EE mDfst - D?s t D?](jf\tf? Q _ + [DfD%s t D%s t Df ](#|<#)£ 

~ [Dfs t D% - s t DtDf]{g?\g%y Q -_ - [D<fs t D? - s t D^ D^]^f\g^f Q+ }. (110) 
From the relation (31) we see that [Y, Hs] = 0. If we define the effective evolution operator by 

V t = {e^ Hs ® l)U t (111) 
which satisfies the quantum stochastic differential equation 

dV t = (-^H s eTn Hs ® dt)U t + (e^ Hs ® l)dZ7 t . (112) 
m 

Explicitly this gives 

dF t = (e^ ffs ®l)[E E {^0'®^^'^)-^'®^^'^)}-^+^^}^]^! 
uef <p,<p'eB 

We see that i/'s is a physical addition to the system Hamiltonian due to the presence of the quantum field. This is 
precisely the Lamb shift. Since Y commutes with Hs it is enough to compute its expectations for eigenstates <fi E B; 

= E E {(4>,DfD%<t>) {9l\9tT Q - + + (<f>,D?D%U) W^TqJ- ( 113 ) 

ugf j,k=l 

This is however equivalent to the (complex) shift one calculates using second order perturbation theory. For example, 
taking the zero temperature for simplicity, one calculates in second order shift |19j 

E E E / dT{4>®^ R ,D^®A^y^- E ^l h Df,®A\gf,)4>®* R ). (115) 

Here we have used the well known identity 

dt eM = 7 1 = nd W ~ f P(-)5 & e M )' ( 116 ) 
-oo i(x-tQ+) x 



where p means that we take the principal part of the integral. Now D"4> is an cigenstatc of Hs with eigenvalue — fko, 
so the summation need only be considered over tu = uj' in (2.6:17) above. Therefore we have 

N(w) 



Y P = EE/ dr(<}>®* Rl Df ®A{g^ H «- h ^/ h Dt®A\gt)<}>®* R ) 
ueF j,k J -°° 

JVH o 

Y,Y.^ D f D ^) / dr^ R: A(g-)AHS-g^ R ). (117) 



wef j,k 

(2) 

Hence Yf> = (</>, F</>). The real and imaginary parts of Y" are therefore the linewidth and energy shift as would normally 
be calculated using second order perturbation theory, this is true in the non-vacuum cases also. 



2.7 Transition Probabilities 

Let Pt(ifj\(j)) denote the probability that the sytem will be measured in state i/j at time t if it initially was in state <f>, then 

p t (V#) = (</> ® Ut{\i>){^\ ® l)C/ t </) ® (118) 

from theorem 3,3a or 3b we have 



d 

Therefore, setting t = 0, 
If ?/; = </> we obtain the relation 
while if (4>,ijj) — then (112) gives 



j t vt{m) = {4>® CZ/iod^^l) ® if/^ ® (119) 

^ft(W)|t=oH&£o(MM)$- (120) 

^W#)|t=o = -<</>,!>), (121) 

^Pt(^)lt=o = <^e(|v)W)^), 

£ Y,{{4>,Df^){i>,Dt<l>){g?\gtY Q+ +(<t>,DW)(Tp,Dfd>}(g?\g%)%_}, (122) 

w6F j',fe=l 



Using the relation dt e lxt = 2-k8{x) we can rewrite this as 



UJ 



where q(k) is the spectral function associated with Q, cf (125). This is our formulation of the Fermi golden rule for 
transitions of the systems state and it corresponds to the usual expressions, cf formulae (1.21. 27a, b) of [2]. 



3 THE WEAK COUPLING LIMIT IN QED. 

As an illustration of our theory we consider the case of quantum electrodynamics. We stress however that the theory 
encompasses a wide range of physical phenomena. For instance exciton models in solid state physics such as phonon 
models or the Frolich [15], [16] polaron model will differ from the following treatment in only minor technical details. The 
electromagnetic field acts as reservoir for our system (S) which we take to consist of a single electron. The electromagnetic 
field can be derived from the potential A given by 



Here we consider two transverse polarizations (a = 1,2) for each mode k. In our notation {e 1 (k), e 2 (k), k = | k| 1 k} 
form a right-handed triad for each k. This ensures that we are working with the radiation gauge V.A = 0. The operators 
%.(k) on the reservoir state space T-Lr satisfy Bose commutation relations 



[a CT (k),4,(k')] =<W5(k-k'). (125) 
The total Hamiltonian for the system and reservoir is 

H = -^|p-eA| 2 + $(r)+^ 
Ira 

= H S + H R + H I +H' I , (126) 
where the unperturbed system Hamiltonian (with potential $(r)) is 

^5 = ^|p| 2 + 1>(r); (127) 
H R = J2 I d3k hc \ k \ 4(k)a CT (k), (128) 

(7=1,2 •* 

-~ E f d3k {4(k)e- lk - r + Mk)e ik - r } G CT (k).p, (129) 

„ — 1 n -J 



Hi 
with 

and 



If we rescale the electronic charge as e Ae we find that 

H ^ H S + H R + XHj + \ 2 H'i. 
In the subsequent analysis we shall drop the term A 2 H' / and consider only 

HW =H s + H r + XHi. (130) 

It has been established rigorously that this does not effect the final result in the weak coupling limit. Now the interaction 

Hi given by (120) has response terms described by the vectors 0?(k) = - — e _tk r Gf (k).p. We assume, as usual, that the 

J nm _ J 

unperturbed system Hamiltonian Hs has a complete orthonormal set of of eigcnstates B. In the case of the Hydrogen 
atom, this means that we consider only the bound states and ignore the effect of the ionized states. In general, Us 
can be decomposed into complementary subspaces generated by the discrete, the absolutely continuous and the singular 
parts of the spectrum of Hs- It is enough to prepare the system in the discrete spectrum subspace to apply our results. 
It is the standard approach in atomic physics to study only the behaviour of bound states anyway, so we are justified in 
this restriction. 

We introduce the test-functions 

.9^(k) = £^(^|e- lk - r p|^)-G CT (k). (131) 

The interaction Hi can be expressed as 

Hi = E f d3k (4(k)5^(k) - a„(k)^(k)} ® T w 

o-=i,2 J 

E E f d3k ® 4(k)5^'(k) - 7> ® a a (k)s^,(k)} 

^ {T^^A^g^-T^^Aig^)}, (132) 



where g^i = gh^, © g\^, € L 2 (R?) © L 2 (R 3 ) = H R and arc the creation/ annihilation operators on T b(Hr) 
®Y B {L 2 {R)); 

AHf 1 ®f)= £ f d 3 kr(k)al(k), 



a=l,2 * 



^J 1 ©/ 2 ) = ]T fd 3 fc/ CT (kK(k). ( 133 ) 

(T = l,2 

Our choice of T-L R above for one particle of the reservoir space is quite natural; namely it consists of wave-functions in 
the momentum representation with two transverse polarizations. The state of the reservoir is in our case determined by 
the covariance operator Q which we shall now specify as that of a thermal field at inverse temperature /3)0, that is 

Q : H R ^ U R : h x © h 2 i-> h\ © h 2 

with Mk) = q{c\k\)h a (k), (134) 

where q(ui) = coth — . With ui^i = [E^i — E^)jh G F,we define 



St**' : U\ ^>H R :h 1 ®h 2 ^h 1 ®h 2 



with 



h a (k) = e i(c|k| -^' )t /i (7 (k). (135) 

In this setup we have allowed for the most general coupling, that is where all the fundamental frequencies F = {tu^ : 
(ft, (ft' £ B} are to be considered. This set is always degenerate in general; however it is important to consider two classes 
of degeneracy arising. The first is the secular class; these are the situations in which degeneracies always arise regardless 
of the spectrum {E<f, : (ft € B} of Hs; they are the pairs ((ft, (ft') an d (V'jV'') which have uj^i = cj^i due to one of the 
following reasons 

1. <ft = <ft = ip = ip , 

2. (j> = (ft , ,ip = 1>'; ((ft^ift), 

3. = ^,0' = /; (eft ^ eft'). (136) 

Any solution to the equation uj^> = oj^m , or equivalently — E^i = E^, — E^i , not of the secular type shall be called an 
extraneous solution. The extraneous solutions are of course dependent on the spectrum of Hs- It is standard procedure 
in the physical literature to assume that such possibilities do not arise however this is a requirement on Hs which cannot 
be made in many important examples. For a particle in a rectangular box, apart from the natural degeneracies arising 
if the ratios of the sides are rational, we also have to consider the fact that the contribution to the energy for the mode 
of vibration rii along the i th -axis is proportional to nf , this means solving the Diophantine equations for the harmonics 

2 2/2 ,2 

n i - m t = n i - m i- 

For the Hydrogen atom we have, apart from the spherical harmonical degeneracies, to consider the integer solutions to 
the Diophantine equations 

1111 

n 2 m 2 n' 2 m' 2 

for the principal atomic numbers. After simple manipulations this leads to the study of the intersection of the algebraic 
projective curve in R 4 ; 

222_ 222_ 222, 2 2 2 _ <-< 
^ ]_ *^ 3 *^ 4. X^X^X^ X-yX^X^ I XyXr^X-^ — u 

with the lattice of positive integers. In the weak coupling limit we obtain the quantum stochastic differential equation 

dU t = E E ^ ® dB$(gw,t) + T^, © dBZbw *)} + Ydt 1 U t 

u£f 4>,4>'eB 

where 



uGf tj>,tj>', 1p,ip'eB 



with Uq = 1. However using the fact that TL,T^' = ((f), tl>)T^'^i etc., we may write Y as 

^=E E {(90>I5^)q+ K + (5^'I^')qI}^- (138) 

In the summation we consider only and tp for which there exists a so that = cj^> , however this is equivalent 
to demanding that u^, = as we always have the identity lo^ = uj^i — uj^a,' . Therefore Y is a linear combination of 
terms with tu^ = and this, in particular, implies that Y commutes with Hs- It is natural in light of this to write 
Y as 

Y = ^ 1/4,4, Ttjnp = ^2 U4>i> (139) 
4>,4>eB <p,4>eB 

where 

V<f>i> = E {(30>I50'v)q+* + (944>'\9m')'q- }■ ( 14 °) 

Now the master equation associated with this problem is from (105) 

= L*(s t ) = -(Fst + s t F f ) + ^ ^ [ T l l p>StT 4n i,>(g 4>4> >\g in i,>)Q + +T Hl s t T^,(g^\g ( j > ^)^_ ]. (141) 

m£f 4>,4>' ,tf>,t(}' eB 

In order to find the general expression for H' s , we see that 



E />E /^^(k) 3 ,vk)^ 



'GS 00 cr=l,2 



But using 5^,(k) = -s£ v (-k) we have 

E f dT E / d 3 fc^(k)^:(k){ 



1° dr V / ^fcg- (k)^7(k){ e -l k ^ (c|k P ± 1 +e ^l ^ c ' k l) _ }e ^« 



,2 



E /° dT E / rf3fc E Me-* k -w>^V k - r ^> 

y B J-oo ff= i2*' 7,7' = 1,2,3 



,2 

xGJ(k)GJ,(k){ e - tc ' k ' t g(c|k ^ + 1 + e ^|k|t g( c l k j) __. } e *W 



«4E /° ^E /** E <*|ei»%^>,e 

lis- T~> " — OO 1 Q ^ A A! lOQ 



fi'EB 00 (7=1,2 " j,j'=l,2,3 

xGJ(k)GJ,(k){e— l k l t g(c|k ^ + 1 +e ^|k|t g( c l k j) zl }. (142) 
We remark that the effect of the response term is as follows; from the commutation relations of r and p we have that 

e - ik - r eifi H V k - r = e" lk r cxp + ^(r))e ik - r = cxp 1 (JE±^!l + W r )), (143) 

ftk.p fr 2 |k| 2 

that is Hs is replacd by Hs H 1 . Now the k dependence in the above expression prevents us from using the 

m 2m 

well-known isotropic identity 

J?l M / SG ' (k)G '' (k) " P^T'"- (144 » 



to calculate y^ as in the dipole approximation. Note that the Lamb shift and the damping coefficients are effected by 

inclusion of the response terms. 

1 i 

The complex shift = —T$ + giving the linewidth and energy shift E'^ for a state <j> £ B can the be written 

as 



>E E (m ^^ + S^'L ^ 

a = l,2 J j',j'=l,2,3 



2ihm 2 ^ J ^ ^^ JL V+(k) -i0+ P-(k)-iO 

where the denominators in the above expression are 



V±(k) =H S + ^V + ± Hc\M 25,. (146) 

m Ira 

This expression has been derived in |20| for the zero temperature case; see however |21j . 
APPENDIX. A: The Traditional Derivation of the Master Equation. 

For the sake of comparison we give the standard arguments used in the derivation of the master equation. This section 
follows closely the development of Louisell [2] . The interaction is taken to be of the form 

H I = Y J D 3 ®F 3 , 

3 

where Dj and Fj act nontrivially on the system and reservoir spaces respectively. We assume that Dj evolves harmonically 
in time under the free evolution with frequency Uj We assume that at time t the system and reservoir are uncoupled, 
that is the density operator p(t) at time zero factors as 

P(0) = P { S) ® P {R) • 

No subscript is required for p^ as we assume that it is invariant under the free-evolution (in particular this is true 
for the choice of a thermal state p^ — e~^ HR ~^ /Tre _/3 '" ff «^ A1 '). We define the reduced system state at time t in the 
interaction dynamics to be the density operator 

The iterated series expansion of st, truncated to second order, is; 

St = so + i j dhTm^iHj), s ® pW] + dh jT dt 2 T rnR [v^(Hj), [v^(Hj), s ® p {R) }}, 

where we have set A = 1. Substituting in for the potential Hj we find 



St = s 

3 



f dti{v^ (Fj)) R e~ iu)3tl [Dj, s ] 
3 Jo 

+ E / * / 1 dt 2 e-^ + ^{[D 3 D k so - D k s Dj](vi°\Fj)vg\F k )) R - [D 3 s a D k - s D k Dj}(v^ (F k )vg> (Fj)) R }, 
j,k Jo Jo 

where (■)r = Tru R [p( R \]. Due to the invariance of the reservoir fields under the free evolution we have that 

(v?\Fj)) R = (Fj) R ; 

(vl \FM 0} (Fk))R = {vE(F J )F k ) R . 
Therefore, if we let t = t\ — t 2 , y = t 2 then we obtain 



st = so 

3 



+ Y,(Fj)R[Dj,s } f e-^dy 
3 J ° 

V / dye-^+^y [ \Te- i ^ T {[DjD k s -D k s D j ](vi \Fj)F k ) R -[DjS D k - So D k Dj}(F k vi°\Fj)) R }. 
„■ t jo Jo 



The approximation procedure is based on the following four steps; step I. One postulates that the contributions coming 
from the sum of all terms higher that second order in the iterated series are negligible, step II. One postulates a hnitc 
autocorrelation time r c such that 

(vW(F j )F k ) R = 0= (FjVpHFJU, 
whenever |t|)t c . Thus for t))r c one may replace the upper limit of the T-integral by +00. This gives 



St = so 



{Fj ) R [Dj j so] /' (ujj ) 



where 



^{[DjDkSa - D k s Dj]wf k - [D jSo D k - s D k D j ]w kj }I t (u) :j + oj k ), 



<k = / e-^{v^{F 3 )F k ) R dr 



u 3,k 

j u 
poo 



and 

I\oj) = [ e- iujy dy. 
Jo 

Step III. For t large with respect to r c one makes the replacement 

z*(w) ^ 

iSiep IV. One postulates that the formulas deduced under the previous assumptions when t is large with respect to r c 
hold also in the limit £ — s> 0; this gives 

ds s t -s 

— = hm 

dt 1 t-co t 

OJj—0 LUj~\-LOk—0 

= [Dj>*o](Fj)R + i\ D 3 D kSQ - D k a Dj]wf ik - [D 3 s D k - sqDuD^w^}. 

3 j,k 

The assumptions leading to this equation have a decidedly ad hoc nature, especially those introduced in steps III and 
IV. The replacement for J* , put in by hand, in step III is precisely what is needed to allow the limit to be taken easily. 
It is instructive to calculate explicitly the master equation in a particular case. We consider as reservoir a free Bose gas 
at inverse temperature 8 and fugacity z = . This can be described by the quasi-free state </?q, on L 2 (M. n ) for example, 
characterized by 

(A\f)A(g)) R = Vq {A\f)A{g)) = (/, ^-±g), 



where 

1 + zerP H R 8 

Q=7 ± ^r=cothf(^-M)- 

1 — ze ' R £ 



We may take H R to be for instance —A. 

We may write the interaction Hj of (31) in the form we are considering with the notations 

Hi = ihJ2i D 3 ® A \9j) ~ h.c} = %, a) ® % a) , 

3 0» 

where we have a summation also over an index a G {0, 1} with the notations 

D m = D h D (jtl) = -D]- F m = A\ 9j ), F {jil) = A( 9j ); 

and consequently 



Wc then have 

wfjoviki) = I dre- i ^{v^{A\g j ))A{g k )) R = I dr <p Q (A\g? 9j )A(g k )) 



- 1 



dr{S^gj, -^— 9 k) = { 9 j\ 9 k)A_ , 
jo L 

Similarly, using the CCR, we find 

/>oo 

«&),(*«,) = / dre^(v^{A{g 3 ))A\g k )) R = ( 9j \ 9k )%- , 
w (ko),(ji) = (ffils*)<3~> 



W (jl),m = (9j\9k)g + , 

while wfj € \ fji^fj = if e = e' as we have (A{f)A{g))n = = (A^(f)A^(g))R. We note that (F^ a ))R = in all cases. 
The master equation then reads 

-^\o = E {[Dj, a D k , a 'So - Dk. a <s D ja \w^ a) {kal) 

-[D jta s D kiOt > - s D Ka >D ha ]w~ ka , ) {ja) } 

= - E {[ D i D t 8 o ~ D k*o D j]wyo),{ki) + [DjDk8 - D k s D]]w+ 1)tm 

-[D}s D k - s D k Dt]W( k0)m - [D jSo Dl - soDlDjjw^ ^}, 
or writing in our notations (and employing the relabeling in terms of the frequency degeneracies as in (31)) 

|io = -EE 



D^o - D%s Df]{g»\gZy Q + 
-[Dfs Q D% s D» k Df}(g»\g%r Q -_ - [Dp D? - s Q D^ D«]Ufyff Q+ } . 

But this is exactly ~^\o = ^o(so), where Lq is given as before. 

APPENDIX B: The Convergence of the Collective Processes to the Noise 
Processes 

The mathematical theory behind the weak coupling limit developed in [1] and subsequent papers is the following. We 
estimate the behaviour as A — > of matrix elements of Uty 3 with respect to collective coherent vectors, that is vectors 

of the form B^fL {f {l) )...B^ x2 {f n) )^ R . This involves substituting u t (Fj) into the series expansion for uty^ and 
examining each of the terms arising. The deep analysis of [1] shows that each term upon normal ordering leads to two 
classes of terms: relevant ones (type I) and negligible ones (type II). The type I terms are exactly those put into normal 
order by commuting time consecutive pairs of reservoir variables, the type II terms account for all others. Following this 
resummation it is shown that the type II terms give vanishing contribution in the limit A — > while the explicit limit for 
the type I terms is calculated; uniform convergence is established, the main technical device used here is one of various 
generalizations of the Pule inequality [3] . 

The independence of the noise processes for different frequencies follows from the next two lemmas. 
Lemma 1: for each u G F let f u , f' u € K u and , T u , S'^ , £ R then 

Jim ( X j S^f^du, A / " S^f'^dv) 

A->0 



JS^/X 2 JS'^/X 2 
E (X[S u ,T u ]i X[S'^,T'j)l 2 (R)(/w|./'uj)w 

(®w€f{X[S„,T u ] ® fu>)\ ® W 'SF (X[S> U ,,T'„,] ® fu,')- 



Proof, the left hand side of (B.l) can be written as a sum over uj,uj' G F of terms 

lim / du du'(f^SZf'^)e^-^ x2 - : 
x ^°Js^ J{S' a ,-u)/\ 2 

by the Riemann-Lesbegue Lemma the terms uj ^ uj' vanish while the uj = uj 1 terms converge by inspection to 

(X[S^,T^]^X[S'^,T'^])L 2 (M)(fuj\f' LU )uj 



Lemma 2: for n G N let f!f> G K u ,x£ } G R, S^ '(T^\ for 1 < k < n and each uj E F , then 



exists uniformly for the x 's and [S, T] 's in a bounded set of R and is equal to 

(*« wi^Mx^x^^ ® /£ 1) ))...^(© weF ( a; L") X[s ,„ )iT („ )] ® /("))) *«). 

For a proof see [10]. 

APPENDIX C: The Quantum Stochastic Differential Equation for U t . 

For convenience we shall consider only one coupling frequency uj so that 

j 

and we have dropped the superscript uj from the operators. Also we shall consider only the Fock (vacuum) case Q = 1. 
We define G\(t) by 

&,G x (t)) = Tr{\ct>m ® S^(/')|* fl ><* fl |B^(/)^ a }, 

where G %s and the second inner product is meant on %s®rB(L 2 (R, we know that the limit limA_ i .o(V'; G\(t)} 
exists and equals 

It is easy to show that this limit has the form (ip,G(t)) where : t i-4 G(f) G is weakly diffcrentiablc. In order to 
obtain a differential equation for G{t) we note that for fixed A one has 



G x (t)) = Tr{|^| ® B^(/')|*«)(*fl|B^(/) I ® " D) ® A^fe))^} 



where 

— A 



o 



Now 

r T/A 2 



1 f ' 



rl>,G x (t)) 



i r (T-t)/X 2 

= El A / (SZf,gj)dv(D]l,,G x (t)) 
j X -A-*)A 2 

where we made the substitution u — t/\ 2 = v. We see that for bounded D this converges as A — > a.e. to 

3 

Next of all the term E\ must be reordered as follows 



3 



I^Tr{|0)(V,| ® B^(f)\* R )(* R \B^(f) [D] ® l)u[f x {l®A{S- /x2 ))} 

3 

\ Tr{\4>)&\ ® B$${f')\*n){* R \B%${f) ® 1)1(1 ® ^/aO). ^]} 



A 

"A + "A- 



In a fashion similar to the calculation of Y\, one easily arrives at 



To evaluate the limit of we note that 



^ \ ft/* 2 ftl f^n-l 

7T,— 1 ^ ^ ^ 

00 . „ t /A 2 ^tj .* n _i 

n=1 %tl Jo Jo Jo 

+H I (t 1 )[(l®A(S t/X 2 9j ),H I (t 2 )...H I (t n )]}. 

It can be shown that in the limit A — > only the commutator involving Hi(t\) contributes. Hence 

A / / / 



hmS 6 A = -lim V(-)"" 1 / rfix / dfe... / dt n 

A^o *_> ^JW io Jo 



XX/AS&: ^3 fe }Tr{|0}{^| ® B^C/Ol^flX^fll^C/) 
but this is the same as 



r(A) 
7 s/A 2 



1™ T2 E / d*(S?/#9j,S? /x *gk)(DlDj1>, G x (t)) = - YtSiWW, D]Z» fc G(t)), 



3,k " u 



where we have used a technical lemma (6.3) of [1] 
We now have 



<V,G(t)> = lim(^G A (i)) = (V,G(0)) + lim / (T x (s) + E x (s))ds 

A^O A-s-0 Jq 



Here we have written for (g\f)q + when Q = 1. 

The quantum stochastic differential equation corresponding to this integral equation is 

dU t = {J^iDj ® dB#( gj) t) - D] g dB%( gj) t)) - J2(9j\9kr-D]D k dt} U t . 

3 j,k 

The generalization to Q > 1 and several coupling frequencies uj is now obvious. 
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